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Recent developments of the high-field magnets and re- 
lated experimental techniques make it possible to explore 
the magnetic phenomena revealed only by a strong mag- 
netic field. A typical example is the dimer compounds 
such as TlCuClg, 1 ' 2 KCuCl 3 3 and BaCuSi 2 6 . 4 These 
compounds consist of dimers each being an antiferromag- 
netically coupled pair of spins. Without a strong mag- 
netic field, spins form singlet pairs and the compounds 
are magnetically inactive. It was suggested 5 ' 6 that when 
a strong magnetic field is applied to such a dimer sys- 
tem, triplet excitations are induced and the system may 
exhibit a phase transition at a finite temperature. This 
phase transition can be most naturally interpreted 7 ' 8 as a 
condensation transition of the excitations that behave as 
bosonic quasiparticles. Indeed, the results of the Hartree- 
Fock (HF) calculation 7 of the diluted Bose gas qualita- 
tively agree with the experimental observations, such as 
the characteristic temperature dependence of the magne- 
tization and the algebraic temperature dependence of the 
critical magnetic field. Results of another recent exper- 
iment 9 also identified the phase transition as the Bose- 
Einstein condensation (BEC). Even a quantitative agree- 
ment was obtained 4 recently for BaCuSi20@ for which 
the temperature dependence of the experimentally mea- 
sured specific heat agreed with the Monte Carlo simula- 
tion results of the effective hard-core boson model, which 
showed a clear 'lambda' peak. Thus, the nature of the 
transition as a condensation transition of the triplet ex- 
citations has been established beyond reasonable doubt. 

However, an unsettling disagreement has been left un- 
solved, in estimates of the critical exponent that charac- 
terizes the temperature dependence of the critical mag- 
netic field. The experiments and the theories both sug- 
gest the algebraic dependence, H C (T) — H c (0) oc T^. 
However, the experimental estimates of range from 1.7 
to 2.0, 2 whereas the theoretical estimate (namely, the HF 
value) is (f> = 1.5. 7 A Monte Carlo simulation 10 was per- 
formed recently on the effective spin model for TICUCI3. 
A temperature dependent effective exponent <j){T) was 
defined such that it characterizes the H c — T curve in a 
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certain finite temperature-range centered at T. The sim- 
ulation result showed that 4>(T) is greater than 1.5 but 
seems to approach 1.5 as the temperature is decreased. 
Moreover, the HF calculation was elaborated recently 11 
based on a more realistic dispersion relation for bosons 
determined by experiments. While it is doomed to yield 
only mean-field values of critical exponents, which obvi- 
ously differ from the correct ones for a finite-temperature 
BEC transition, it still provides a fairly good approxima- 
tion for the value of the critical temperature. Therefore, 
these recent simulation and theoretical studies may be 
suggesting that the failure of the HF approximation in 
the finite-temperature critical phenomena does not nec- 
essarily mean the failure in the quantum critical phe- 
nomena. In particular, the mean-field value <j) — 1-5 may 
be the correct and exact value near the zero tempera- 
ture. This statement, however, has not been made con- 
clusive because of the uncontrollability of the approxima- 
tion and/or the various limitations (such as the system 
size and the temperature range) in the Monte Carlo sim- 
ulation. 

In this letter, we present a simple theoretical argu- 
ment, from which one can conclude that the mean- field 
values of the critical exponents are exact for the quantum 
critical point (QCP), i.e., T = 0, H = H c (0), and there- 
fore the inaccuracy of the HF approximation is not the 
reason for the disagreement. To demonstrate the validity 
of the argument, we also show some results of a quan- 
tum Monte Carlo simulation on the three-dimensional 
XY model with a magnetic field. 

We start with the model in which each spin (e.g., each 
Cu spin in TICUCI3) is treated explicitly: 

H = ^Js ll - S i2 + J'ia,jl3 8 ia ' «j/9 ~ ^ Hs z ia . 

i (ij) a,/3 ia 

Here, J > is the antiferromagnetic coupling constant 
that binds a pair of spins, i the index of the dimer and 
Sn and Si2 the two spins that form the pair i. The lat- 
tice structure and the magnitude of the coupling con- 
stants depend on the compound. In the case of TICUCI3, 1 
for example, the lattice is composed of double chains of 
spins (i.e., two- leg ladders) with relatively weak interlad- 
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der couplings. The dominant couplings are theose corre- 
sponding to rungs of the ladder, which are represented 
by J above. Below, we assume that the interdimer cou- 
plings J- a are much smaller (in magnitude) than J, 
and that the lattice is three-dimensional. Real materi- 
als, such as TICUCI3, are often strongly anisotropic. The 
spatial anisotropy is, however, irrelevant in the present 
letter since we focus on the critical properties, in par- 
ticular, the critical properties near T = 0, where the 
three-dimensionality manifests itself. 

Following Tachiki and Yamada, 5 ' 6 we neglect the two 
of the triplet states that are not favored by the magnetic 
field. The remaining two states, the triplet state that 
is favored by the magnetic field and the singlet state, 
are regarded as the up and down states, respectively, of 
an effective S — 1/2 spin. The resulting Hamiltonian in 
terms of the effective spin Si is 5,6 

h = J2 h (sfsj + svs? + J2 HiS?, (i) 

to) V J i 

where and Hi can be expressed as simple linear com- 
binations of J, H and J[ a jp- We here assume that the ef- 
fective Hamiltonian is translationally invariant and non- 
frustrated. The experimental estimations of the lattice 
structure and the coupling constants suggest that these 
assumptions are valid, at least approximately, for several 
compounds such as TICUCI3 12 and BaCuSi206- 13 

Thus, we have arrived at the uniform 5 = 1/2 XX Z 
model in three dimensions with a spatial anisotropy and 
an easy-plane spin anisotropy. The common belief is that 
the spatial anisotropy does not affect the critical prop- 
erties of the model both at the finite temperature tran- 
sition and the QCP. In addition, the S z -S z term can be 
dropped without affecting the critical properties, since it 
does not change the type of spin anisotropy. Therefore, 
we can further simplify the model without changing the 
critical properties: 



(2) 



to) 



i.e., the XY model Hamiltonian. The sign of the cou- 
pling constant J can be switched by changing the spin 
representation bases. Therefore, we assume J > and 
H > without loss of generality in what follows. 

If we consider the effective Hamiltonian as the sum of 
pair Hamiltonians H^ , it is straightforward to show that 
the state in which the two spins Si and Sj are aligned 
with the magnetic field is the unique ground state of 
Hij when H > dJ. Therefore, the state with all spins 
pointed up is the unique ground state of the entire system 
for H > dJ. On the other hand, it can also be shown 
easily that mixing with one-magnon states decreases the 
energy of the all-up state when H < dJ. Therefore, for 
the present system, the quantum critical point is located 
exactly at H — dJ. 

While the Hamiltonian (2) is the one that we simulate 
with the quantum Monte Carlo method, we can further 
reform it by using the bosonic representation = b\, 
S^ = hi, and Sf = fn — \ with the constraint 
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Fig. 1. Field dependence of magnetization p = i — (S z ). The 
three curves correspond to different system sizes, while the 
'scaled' temperature L 2 / (3 is fixed to be 8. The inset is the close- 
up view of the quantum critical region. 



Then we obtain 



to) 



(b\bj + bth) 



rii 



h^hi(hi - 1), 



(4) 

with new parameters t and /j, corresponding to the old 
ones as t ~ J/2 and /i ~ H. The A-term imposes the 
condition (3), approximately. 

Generally, the replacement of the original 'hard' con- 
straint by a 'soft' constraint may affect the critical prop- 
erties around the QCP. In the present case, however, we 
consider that this modification does not change the es- 
sential properties because the typical occupation number 
is zero or much smaller than unity near the QCP even in 
the soft-core case, while the difference between the hard 
and soft constraints manifests itself only for the states of 
which rii > 2. 

In the following few paragraphs, we present the re- 
sults obtained earlier 14 ' 15 for the Bose gas that can be 
found also in textbooks. 16 We reproduce them here since 
it is convenient to have explicit formulas in the form di- 
rectly comparable with the present Monte Carlo results. 
(In particular, explicit formulas for the system-size de- 
pendence are hardly found in the literature.) First, the 
continuous field theory for model (4) is characterized by 
the action 



hi 



dr 



— hReip — r \ip\ 2 + u \ip\ 



(5) 



The parameter r depends on [i as r oc \i — fi c near the 
QCP, and the symmetry-breaking field h, which does not 
correspond to the original magnetic field, is included for 
a technical purpose. From a simple dimensional analysis, 
the effect of the scale transformation up to the scale b 
turns out to be the following: 



fi- 
ll- 



(3 = fib- 2 , L 
h = hb 2+ i, r 



h t = b\bi = 0, 1. 



L = Lb~\ i> i> = #2 
> f = rb 2 , u — > u = ub 2 ~ d , 



(3) 
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Fig. 2. Size dependence of susceptibility at quantum critical 
point, H = 3J. Three values of the scaled inverse temperature, 
(3/L 2 = 8, 16 and 32, are examined. The solid lines indicate the 
slope of x <x L 2 - 5 . 



Fig. 3. Finite size scaling plot of susceptibility at j3 = 4J. The 
critical exponents used here are taken from the results of the 
classical XY model. 17 



where L is the system size. Hence, the upper critical di- 
mension d c is equal to 2. 14 For three dimensions, there- 
fore, we should expect a mean-field-type scaling. 

Because of the dangerous irrelevant parameter u, the 
usual procedure of the finite size scaling, i.e., rescaling 
the coordinate and the abscissa by size-dependent fac- 
tors, does not work. We first define $(h,r,u, (3, L) as 
the singular part of — log Z(h, r, u, (3, L). We assume that 
$ has the property &(h,r,u, /?, L) — &(h,r,u, /3,L) for 
an arbitrary b. With this expression, we first consider 
the particle density p, which corresponds, in the original 
model, to the deviation of the magnetization (parallel to 
the magnetic field) from its saturation value. The parti- 
cle density is related to the magnetization perpendicular 
to the field, which we denote as m, by p ~ m 2 ~ \ip\ 2 . 
Then we have 



m(r, u) 



lim lim lim L (3 

ft— >0 L^oo /3^oo Oh 



b 2 m(r, u), 



where m(r, u) is the magnetization of the renormalized 
system. When the system is renormalized up to the point 
where the correlation length is 0(1), we can neglect fluc- 
tuations and apply the mean-field-type calculation with- 
out a serious error in estimating rh, which yields 

rh(f,u) <• 

v a y u 

for sufficiently large b. Thus, we arrive at 14 

pocr oc \H-H C \ (T = 0). (6) 

Similarly, we can obtain the scaling form of the sus- 
ceptibility (of the spin components parpendicular to the 
magnetic field) at r — 0: 




,(3,L) = lim lim / dx f dr (i/)*(x,tU(Q,0)) 

J a 



cx b 2 x(uJ,L). 
For sufficiently large 6, we have 
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Fig. 4. Critical field as a function of temperature. The solid curve 
corresponds to <j> = 1.5, whereas the dashed curve <fi = 2.0. 



where L is assumed to be larger than b. Therefore, the 
asymptotic system size dependence of the critical suscep- 
tibility for a fixed value of the scaled inverse temperature 
$ = P/L 2 is 



X {u,(3'L 2 ,L)KL^i. 



(7) 



Finally, we consider the temperature dependence of 
the critical magnetic field. When d > 2, we can make u 
small by choosing large b so that the perturbation theory 
in u is nearly exact for the renormalized system. The 
perturbation theory indicates that $ for the renormalized 
system has a singularity when the condition 16 f = f c = 
Au/(3 d l 2 is satisfied with A being a numerical constant. 
This means that in terms of the bare coupling constants, 
the singularity occurs when 14 

H c (T)-H c (0)ocT d ^ 2 . 

We now present the results of our Monte Carlo sim- 
ulation for spin Hamiltonian (2) on a simple cubic lat- 
tice and compare them with the theoretical predictions. 
The simulation method employed here is based on the 
directed-loop algorithm 18 for which a review can be 
found in a recent article. 19 The system sizes that we ex- 
plored range from L = 4 to L = 32 and the temperatures 
from j3 = 0.5 to j3 = 32 (for L = 32). 
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Fig. 5. Logarithmic plot of critical field against temperature. The 
two straight lines correspond to <j> = 1.5 (solid) and = 2.0 
(dashed). In the intermediate temperature region, </> = 2.0 seems 
to fit the data well, while in the low-temperature region, the 
correct slope cf> = 1.5 yields a better fitting. 



Figure 1 shows the magnetic field dependence of the 
deviation of the longitudinal magnetization from its sat- 
uration value. In the bosonic language, this is the par- 
ticle density which is a function of chemical potential. 
Results are shown only for three system sizes L = 4, 8 
and 16 with fixed 'scaled' temperature, L 2 //3 — 8. As 
can be seen in the figure, the system size dependence re- 
mains only in the vicinity of the quantum critical point. 
The inset shows the magnified view of the critical region. 
The linear field dependence (6) can be clearly observed 
in this inset. 

Next, we show the susceptibility of the perpendic- 
ular components of spins. This is defined as \ = 
J2 r Jo dT(S x (r,T)S x (0,0)) for model (2). In Fig. 2, the 
system size dependence of the susceptibility at the criti- 
cal point H = 3Jis shown. Three different values of the 
scaled temperature L 2 / (3 are chosen. For each value of 
the scaled temperature, the asymptotic size dependence 
is well fitted by the straight line with the slope 2.5, in 
agreement with (7). Thus the universality class of the 
QCP is identified as the mean field type. 

In order to obtain the phase diagram in the H — T 
plane, we carry out the finite-size scaling of the suscep- 
tibility with a fixed magnetic field or a fixed tempera- 
ture. As an example, a finite-size-scaling plot is shown 
for (3 J = 4 in Fig. 3. The data can be scaled nicely with 
the critical exponents estimated for the classical XY 
model, namely, v = 0.67155(27) and r\ = 0.0380(4). 17 
The value of the critical field is the only free parame- 
ter determined by the present data, and the best plot 
is obtained for H C (T = 4J) = 2.862. In this manner, 
we estimated the critical field (temperature) at various 
temperatures (fields). 

In Fig. 4, two curves are plotted for comparison. The 
solid curve represents the mean-field critical exponent, 
4> = 1.5, whereas the dashed curve represents the pre- 
vious estimate 2 based on an experiment on TlCuCbj, 
4> = 2.0. At a first glance, it seems that the curve with 



4> = 2.0 fits the data better than that with (f> = 1.5. 
However, when the logarithmic scale is used, as we do 
in Fig. 5, it is rather clear that (f> — 2.0 explains only a 
transient behavior, and the correct asymptotic value of 
the exponent is = 1.5. 

To summarize, we have argued that the mean-field ex- 
ponents should be correct and exact for the quantum 
critical point of the XY model in three dimensions, and 
our numerical simulations demonstrated that it is indeed 
the case beyond reasonable doubt. Therefore, we can at 
least exclude the inaccuracy of the Hartree-Fock approx- 
imation as the source of the discrepancy between the 
theoretical and experimental results. We have also found 
that there is a relatively large intermediate temperature 
region where a transient behavior can be well described 
by the effective exponent 0^2. While there could be 
some other possibilities that may truly change the criti- 
cal behavior from the mean field type to something else, 
such as a presence of Dzyaloshinsky-Moria interactions, 
the present results suggest that it is possible to explain 
the discrepancy between the theoretical and experimen- 
tal results only by taking into account the transient be- 
havior and that the quantum critical behavior driven by 
the magnetic field can be correctly described by the XY 
model (or the diluted Bose gas fixed point) which yields 
the mean- field exponents. 
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